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Abstract. We analyze entanglement generation between a pair of neutral two level 
atoms that are initially excited in a common electromagnetic vacuum. The nonlocal 
correlations that appear due to the interaction with the field can become entanglement 
when the field state is known. We distinguish two different situations: in the first, the 
field remains in the vacuum state, and second two photons are present in the final state. 
In both cases we study the dependence of the entanglement with time and interatomic 
distance, at ranges related with locality issues. 

PACS numbers: 03.67.Bg, 03.65.Ud, 42.50. Ct 



As is well known, an atom B cannot become excited by a physical signal coming 
from other atom A placed at a distance r at a time t < tr = r /c. This does not prevent 
the emergence of correlations between these two atoms aX t < t^. In particular, when 
the final state of the field (and/or atom A) is fixed, atom B can become excited even 
when both atoms were initially (at t = 0) uncorrelated and the time elapsed was smaller 
than tr- The specification of this final state is highly nonlocal and a finite excitation 
probability Pb (t) for atom B at times t earlier than tr does not imply causality violation 
but is an artifact of that nonlocality. As shown in [1], causality is fully restored after 
summing over all the accessible final states for atom A and the electromagnetic field. 
Now only the state of B is accounted for; this is a local property, no wonder that 
Psit) = for t < tr in this case. In the language of Quantum Information, the system 
is nonsignaling [2]. It is well known that the mere existence of correlations is not 
enough in order to send information between the parties: classical communication is 
always needed to transform entanglement into information. 

Due to the existence of these nonlocal correlations, the reduced state of the atoms 
at a time t is a mixed classically correlated state, even for c t < r |3] . The generation of a 
correlated state from an uncorrelated one only by local interactions could be blamed on 
the fact that the field vacuum is an entangled state between spacelike separated regions 
[H El E] or on the finiteness of the Feynman propagator for ct < r [|7j . There are at least 
two ways to obtain entanglement from these correlations: introducing a suitable time 
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dependent coupling between the field and the atoms (like in [HI E] with a scalar field 
and two classical detectors) or fixing the field state |3l El [10]. Here we will show new 
examples of the latter case: we will start with both atoms excited in the electromagnetic 
vacuum, and after a interaction time t we will consider the following field states: 

i) no photons are present, 

ii) two photons are detected. 

We shall study the dependence of the entanglement with t and r, both for r > ct 
and r < ct. 

We assume that the wavelengths relevant in the interaction with the atoms, and 
the separation between them, are much longer than the atomic dimensions. The dipole 
approximation, appropriate to these conditions, permits the splitting of the system 
Hamiltonian into two parts H = Hq + Hj that are separately gauge invariant. The first 
part is the Hamiltonian in the absence of interactions other than the potentials that 
keep A and B stable, Hq = Ha + Hb + -f^field- '^^^ second contains all the interaction 
of the atoms with the field 

Hi = -- d„(x„,t)D(x„,t), (1) 

^0 n=A,B 

where D is the electric displacement field, and d„ = I]j e / d?^i ( | (x^ — x„) | G ) 
is the electric dipole moment of atom n, that we will take here as real and of equal 
magnitude for both atoms (d = = de), | E ) and | G ) being the excited and ground 
states of the atoms, respectively. 

In what follows we choose a system given initially by the product state, \ 4')o = 
\EE)-\0)m which atoms A and B are in the excited state | E ) and the field in the 
vacuum state | ). The system then evolves under the effect of the interaction during a 
lapse of time t into a state: 

|^)* = e-^/o'^*'^^(*')/^|V^)o (2) 
that, to order a, can be given in the interaction picture as 

|atomi,atom2, field) t = {{1 + a) \ E E) + b\G G)) \ 0) 

+ {UA \GE)+UB\EG))\l) + {f\EE)+g\GG))\2) (3) 

where 

a = liOmS^SA + 5+5^)10), b = (0|r(5^5^)|0) 

UA = {l\SX\0),UB = {1\Sb\0) (4) 

/ = \{2\T{SiSj + S^SbM, g = {2\T{SBSjm, 

being S = — f /o dt' Hi{t') = 5+ + S^, T being the time ordering operator and 
\n), n = 0, 1, 2 is a shorthand for the state of n photons with definite momenta 
and polarizations, i.e. | 1) = |k, eA), etc. The superscript signs in 5^, to be defined 
in [5] below, are associated to the energy difference between the initial and final atomic 
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states of each transition. Among all the terms that contribute to the final state only b 
corresponds to interaction between both atoms. This would change at higher orders in a. 
Here, a describes intra-atomic radiative corrections, ua and ub single photon emission 
by one atom, and g by both atoms, while / corresponds to two photon emission by 
a single atom. In Quantum Optics, virtual terms like b and / which do not conserve 
energy and appear only at very short times, are usually neglected by the introduction 
of a rotating wave approximation (RWA). But here we are interested in the short time 
behavior, and therefore all the terms must be included, as in [U IHl |12] . Actually, only 
when all these virtual effects are considered, it can be said properly that the probability 
of excitation of atom B is completely independent of atom A when r > ct im 112] . 
Finally, in the dipole approximation the actions hS^ in (jl]) reduce to 

5± = -M* rft'e±^™'dE(x,t') (5) 
n Jo 

where Q = uje — ojg is the transition frequency, and we are neglecting atomic recoil. 
([5]) depends on the atomic properties O and d, and on the interaction time t. In our 
calculations we will take (f2|d|/ec) = 5 • 10^^, which is of the same order as the Is — >■ 
2p transition in the hydrogen atom, consider f2t > 1, and analyze the cases (r/ct) ~ 1 
near the time t = t^ where one atom could begin to receive signals from the other. 

Given a definite field state | n ) the pair of atoms is in a projected pure two 
qubits state as shown in ([3]). We will denote these states by \A,B,n), = 
\A, B, n) (A, 5, n \ and will compute the concurrence C^"-* [13j : 

C(") = max{0, ^fXi-Y. V^} (6) 

where A, the largest of the eigenvalues (j = 1, 4) of [(a^ ® (^y)PAB i'^y ® '^y)]pAB- 
We begin with the case n = 0, where the field is in the vacuum state and, following 
([3]), the atoms are in the projected pure state 1^450) = ((1 + a) \ E E ) + b \ G G )) / Cq, 



where Cq = y |1 + op + |&P is the normalization, giving a concurrence 

= 2\b\\l + a\/cl. (7) 

The computation of a and b can be performed following the lines given in [3], where 
they were computed for the case of a initial atomic state \ E G). We will consider that 
the dipoles are parallel along the z axis, while the atoms remain along the y axis. Under 
that conditions, using the dimensionless variables x = r/ct and z = Qr/c: 

AiKz^ Zr, 



3 X 



n |1 I + 2 Z tt). 



^(-V% + V.V,)/, (8) 

o2 ^2 



with K = a\dY' / (e r^) and 1 = 1+ + /-, where: 



Ei{tz{l±-)) - e''^^^^^ Et{-tz{l±-))] (9) 
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for X > 1 with the additional term — 2 vri e*^*^^"^/"^'^ for x < 1. We use the conventions 
and tables of |14j . 

We show in Fig. 1 the concurrence C^^-* (171) for x near 1 for given values of z. Like 
the case where \EG) is the initial atomic state, C^^-* jumps at a: = 1 and has different 
behaviors at both sides. 
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Figure 1. Concurrence C'^°-' of the atomic state in the e.m. vacuum as a function 
of a; = {r/ct) for z ~ {flr/c) = 5 (solid hue), 10 (dashed line) and 15 (dotted line). 
The height of the peak is C'^°-' = 1. x ^ (t ^ oo) is the region usually considered in 
Quantum Optics. 



In Fig. 2 the concurrence is sketched as a function of z for given values ofQt = z/x. 
The tiny values of the concurrence for the region z > Qt (which corresponds to x > 1), 
diminish as t grows and will eventually vanish, since 6 is a non-RWA term. 
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Figure 2. Concurrence C^°-' of the atomic state in the e.m. vacuum pj^^ as a function 
of z = (ilr/c) ioi z/x ~ ilt = 6 (solid line), 9 (dashed line) and 12 (dotted line). 
X > 1 amounts to z > 51 1 in each case. 



The case n = 1 analyzed in p], shows atom-atom correlations due to the uncertainty 
of the photon source. Here, we will focus on the two photon case. The final atomic state 
\ AB2) = {f \E G) + g \ G E))/c2, with C2 = ^| / P + I fl' is in the same subspace 
as for n = 0. The normalization C2 is 0{a) like the expectation values /, g, so that all 
the coefficients in p^"^^ may be large. Therefore, although the probability of attaining 
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this state is small, the correlations are not. The concurrence is 

C(2)=2|//|/c^. (10) 

We find that: 

f = e{ti-t2)ivAitl)uAit2) + UA{tl)v'j,{t2) 

+ VB (ti) u'b (t2) + UB (ti) v'b (t2) ) (11) 

g = ubu'j^ + ua u'b 

with va = ( 1 I 5 I ) and vb = (l|iS^|0). The primes account for the two 
single photons, i.e. |2) = |keA, k'cA'). The quantities | ma P = I p = | m P, 
I f A P = I f B P = I ^ = uav*^ = ubv% uv* = u\v\ = ubv*^, ubu\ and 
Va v*b have been computed in ^ . 

In Fig. 3 we show C^"^^ in front of x for given values of z. When x ^ (t oo, 
i.e. the Quantum Optics regime), / vanishes and the final atomic state would be the 
separable state \GG), with zero concurrence. Entanglement is sizeable for a; > 1, and 
could be maximized if a particular two photon state was detected |10j . 
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Figure 3. Concurrence C'^' of the atomic state with n = 2 photons p^^B front of 
X — r/{ct) for z = flr/c = 5 (sohd hne), 10 (dashed line) and 15 (dotted line). 

In Fig. 4, C^^-* is sketched as a function of z for given values of Vlt = z/x. Again, 
the concurrence for the region z > Qt [x > 1), diminish as t grows and will eventually 
vanish, since it is due to /, which is a non-RWA term. Interestingly, as we noted in [9] 
for the single photon emission, x = 1 is a singular point that divides the spacetime into 
two different regions. This occurs, as in [9], even if in these cases t is not the propagation 
time of any physical signal between the atoms. This effect comes from the appearance 
of effective interaction terms like /, that would be missing if we could discriminate the 
source of emission of each photon. 

As a summary, we have computed the concurrence generated between a pair of 
neutral two level atoms by the interaction with the electromagnetic field, in order to 
study the role of locality in the growing up of quantum correlations. We have considered 
the initial state with both atoms excited and the field in the vacuum, and two different 
situations after a time t: the field remains in the vacuum and two photons have been 
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Figure 4. Concurrence C'-^^ of the atomic state with n = 2 photons in front of 
z = Qr/c for fit = z/x — 6 (sohd hne), 9 (dashed line) and 12 (dotted line). 

emitted. In the former case, entanglement is generated by the interaction amphtude 
b. For r < ct photon exchange make correlations grow during a short time, before 
vanish eventually due to the non-RWA nature of b. For r > ct although b is nonzero 
its amplitude is negligible compared with r < ct. The situation is different when two 
photons are present. There is a relevant probability that both photons come from the 
same atom, which in our scheme is represented by the non-RWA amplitude g. As long 
as we cannot distinguish the source of emission of each photon, concurrence is sizeable 
for r > ct and r < ct, and will eventually vanish. Although there are no interaction 
terms between the atoms, the indistinguishability construes effective interaction terms 
and r = ct plays the role of a frontier between two different spacetime regions. 
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